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ABSTRACT How many distinct molecular paths 
lead to the same phenotype? One approach to this ques- 
tion has been to examine the genetic basis of conver- 
gent traits, which likely evolved repeatedly under a 
shared selective pressure. We investigated the conver- 
gent phenotype of blue iris pigmentation, which has 
arisen independently in four primate lineages: humans, 
blue-eyed black lemurs, Japanese macaques, and spider 
monkeys. Characterizing the phenotype across these 
species, we found that the variation within the blue- 
eyed subsets of each species occupies strongly overlap- 
ping regions of CIE L*a*b* color space. Yet whereas 
Japanese macaques and humans display continuous 



variation, the phenotypes of blue-eyed black lemurs 
and their sister species (whose irises are brown) occupy 
more clustered subspaces. Variation in an enhancer of 
OCA2 is primarily responsible for the phenotypic differ- 
ence between humans with blue and brown irises. In 
the orthologous region, we found no variant that distin- 
guishes the two lemur species or associates with quan- 
titative phenotypic variation in Japanese macaques. 
Given the high similarity between the blue iris pheno- 
types in these species and that in humans, this finding 
implies that evolution has used different molecular 
paths to reach the same end. Am J Phys Anthropol 

151:398-407, 2013. © 2013 Wiley Periodicals, Inc. 



In natural populations, the number of possible genetic 
changes leading to a given trait is determined by the 
combination of two factors: the number of loci at which 
mutations influence the phenotype (the mutational tar- 
get size) and the proportion of these mutations that do 
not have prohibitively deleterious consequences on other 
phenotypes (the extent of pleiotropy; Stern and Orgo- 
gozo, 2008; Stern, 2011). Examples of convergent evolu- 
tion, the independent acquisition of the same trait in 
different populations or species, may be highly informa- 
tive with regard to questions of evolutionary predictabil- 
ity and constraint because these likely represent 
multiple solutions to similar selective pressures (see, 
e.g., Conway Morris, 2003; Christin et al., 2010; Losos, 
2011). 

To date, cases of convergence have been found to occur 
through independent coding changes in the same gene 
(e.g., Mundy, 2005; Christin et al., 2008; Gross et al., 
2009; Kingsley et al., 2009); through changes in distinct 
non-coding elements regulating the same gene (e.g., 
Miller et al., 2007); through distinct changes in the 
same regulatory element (e.g., Tishkoff et al., 2007); or 
even through changes in the same protein domains (e.g., 
Aminetzach et al., 2009; Feldman et al., 2012) or at the 
same amino acids (e.g., Li et al., 2010; Liu et al., 2010b; 
Shen et al., 2012). In other cases, changes in distinct loci 
result in the evolution of the same trait (e.g., Chen 
et al., 1997; Hoekstra et al., 2006; Borowsky, 2008; 
Steiner et al., 2009). Most studies have found that con- 
vergent phenotypes tend to evolve through changes 
within the same gene (Conte et al., 2012; Martin and 



Orgogozo, in press) or within a particular pathway (e.g., 
Chan et al., 2012; Tenaillon et al., 2012); in the absence 
of a reporting bias, these findings indicate that the num- 
ber of permissible paths to a trait is often limited. 

Many of these reported cases of phenotypic conver- 
gence involve pigmentation (see Hubbard et al., 2010), 
notably in the color of skin (Miller et al., 2007; Gross 
et al., 2009), fur (Hoekstra et al., 2006; Kingsley et al., 
2009; Steiner et al., 2009), or plumage (Mundy, 2005). In 
primates, blue iris pigmentation has been documented 
in four different lineages. In humans (Homo sapiens), 
the prevalence of the trait increases with latitude from 
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24 to 55% in European populations (Zanetti et al., 1996; 
Laeng et al., 2007); it has also been observed in popula- 
tions with European admixture (Frudakis et al., 2007). 
In Japanese macaques {Macaca fuscata), which are 
endemic to the broad-leaved deciduous and evergreen 
forests of Japan (Oi, 2002; Abe et al., 2005), blue or in- 
termediate iris color has been observed at frequencies of 
12 to 19% in colonies on Shodoshima and Kyushu 
islands (Yamagiwa, 1979; Zhang and Watanabe, 2007). 
Blue irises have also been reported in the brown spider 
monkey (Ateles hybridus, formerly Ateles belzebuth 
hybridus; Hernandez-Camacho and Cooper, 1976; Konst- 
ant et al., 1985) and closely related Colombian black spi- 
der monkey (Ateles geoffroyi or fusciceps, subspecies 
rufiventris; Hernandez-Camacho and Cooper, 1976; 
Defler et al., 2004), which inhabit evergreen, semi-decid- 
uous, and montane seasonal forests in Colombia, Pan- 
ama, and Venezuela (Mondolfi and Eisenberg, 1979; 
Defler et al., 2004). Although the prevalence of this trait 
in these species has not been extensively documented, 
approximately 13% of a population of brown spider mon- 
keys in north central Colombia has blue irises (R. Rim- 
bach and A. Link, personal communication). In contrast 
to these other species, all blue-eyed black lemurs (Eule- 
mur flavifrons, formerly Eulemur macaco flavifrons) 
have blue irises, whereas in the most closely related spe- 
cies, the black lemur (Eulemur macaco, formerly Eule- 
mur macaco macaco), all individuals have brown irises 
(Mittermeier et al., 2006). These sister species inhabit 
primary and secondary tropical sub-humid forests in a 
narrow range in northwestern Madagascar (Rabarivola 
et al., 1991; Andrianjakarivelo, 2004; Randriatahina and 
Rabarivola, 2004) and hybridize across part of this range 
(Rabarivola et al., 1991). Given that almost all other pri- 
mates have brown or yellow irises (Kobayashi and Koh- 
shima, 2001), blue iris pigmentation can be inferred to 
be a derived trait that has arisen independently on these 
lineages. 

Beyond these broad strokes, the extent of iris pigmen- 
tation variation within and between species has not 
been characterized. In particular, it is not yet clear 
whether the reported blue irises represent the same phe- 
notype in all species. Demonstrating such similarity is 
important, as if the derived phenotypes were different, 
we would anticipate the involvement of distinct genetic 
loci. Thus, mapping the traits would not provide addi- 
tional information with regard to evolutionary con- 
straint. Quantitative measures of iris pigmentation 
variation, recently demonstrated to be associated with 
genetic variation in humans (Liu et al., 2010a; Edwards 
et al., 2012), provide a way of evaluating phenotypic 
similarity across species. 

In humans, the phenotype of blue iris pigmentation 
has well understood molecular and genetic underpin- 
nings. Observable human "eye color" depends upon the 
amount and type of melanin within the outermost layers 
of the iris (Prota et al., 1998; Wielgus and Sarna, 2005). 
Whether human irises are blue or brown can be pre- 
dicted with >80% accuracy by the genotype at a single 
nucleotide polymorphism (SNP), rsl2913832, in a regu- 
latory element influencing the expression of OCA2 
(Eiberg et al., 2008; Sturm et al., 2008; Liu et al., 2009; 
Visser et al., 2012). The 0CA2 gene itself encodes a 
transport protein expressed in the iris that influences 
the development of melanocytes (Rinchik et al., 1993; 
Brilliant, 2001). In addition to the main effect SNP, var- 
iants within OCA2 and at least 12 other genes have 



been more weakly associated with human iris color vari- 
ation; these likely act as modifiers of iris color with 
smaller effect size (Kanetsky et al., 2002; Frudakis 
et al., 2003; Graf et al., 2005; Duffy et al., 2007; Kayser 
et al., 2008; Sturm et al., 2008; Liu et al., 2009). Age 
strongly influences human iris color variation, and sex 
has a significant but weak effect (Liu et al., 2010a). 

Our understanding of the genetic basis of blue iris pig- 
mentation in humans can inform hypotheses regarding 
the genetic basis in other primates. In humans, several 
coding sequence mutations and deletions of coding 
regions in OCA2 cause oculocutaneous albinism, or 
severely reduced pigmentation of the hair, skin, and 
irises (Manga and Orlow, 1999; Oetting and King, 1999), 
whereas regulatory SNP rs 129 13832 has a more moder- 
ate influence on hair and skin phenotypes (Sulem et al., 
2007; Branicki et al., 2009). If dramatic pigmentation 
reduction in skin or fur were deleterious in other prima- 
tes, we may expect blue irises to evolve via regulatory 
mutations in these species as well. Notably, SNP 
rs 129 13832 lies within a stretch of 410 bp demonstrating 
strong conservation among mammals; thus, this region 
is a strong candidate for variants with the capacity to 
influence OCA2 expression. Bradley et al. (2009) 
sequenced 166 bp of the orthologous region and found no 
fixed differences between the two lemur species. They 
identified one variant that was polymorphic within the 
blue-eyed black lemurs but not the black lemurs 
sequenced, a pattern that could be suggestive of this 
variant or a linked site contributing to pigmentation 
differences. 

Several lines of evidence support the idea that blue 
irises have fitness effects in humans. The rarity of obser- 
vation of the phenotype in nature, in combination with 
observed deleterious pleiotropic consequences in humans 
(Eagle Jr, 1994; Imesch et al., 1997; Metallinos et al., 
1998; Santschi et al., 1998; Yang et al., 1998; Smith 
et al., 2000; Pingault et al., 2010) and domesticated ani- 
mals (Juraschko et al., 2003; Geigy et al., 2007; Haus- 
wirth et al., 2012), suggests that some benefit would be 
required to overcome purifying selection against these 
pleiotropic effects. Population genetic evidence also sup- 
ports recent positive selection at the locus responsible 
for the majority of the phenotypic variation in humans 
(Voight et al., 2006; Donnelly et al., 2012; Yang et al., 
2012). Although the mechanism driving this population 
genetic signature is unclear, direct selection for blue 
irises to improve short wavelength perception (Born- 
stein, 1973; Laeng et al., 2007), selection for light skin 
pigmentation to increase vitamin D absorption (see, e.g., 
Jablonski and Chaplin, 2000; Parra, 2007), and sexual 
selection (Darwin, 1871; Diamond, 1992; Aoki, 2002; 
Frost, 2006; Laeng et al., 2007) may have played a part. 

Evidence for an adaptive value of blue irises in non- 
human primates is more limited; however, pleiotropy 
and sexual selection may also have played a role in 
these species. In Japanese macaques, the presence of 
blue irises is moderately correlated with reduced pig- 
ment in fur (Yamagiwa, 1979), so selection on lighter fur 
could also have increased the frequency of blue irises. Of 
note, Japanese macaques have the northern-most range 
of any non-human primate (30°21' to 41°08'; Oi, 2002), 
so selection either for increased absorption of UVB or for 
short wavelength perception could have occurred in this 
species. In contrast, neither blue-eyed black lemurs 
(which are more lightly pigmented than black lemurs 
among females; Mittermeier et al., 2008) nor either 
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species of spider monkey with blue irises live at extreme 
latitudes. Several adaptive hypotheses involving blue 
irises, including sexual selection and species recognition, 
have been proposed for the blue-eyed black lemur (Brad- 
ley et al., 2009). Alternatively, blue irises may have risen 
to high frequency by chance (i.e., via genetic drift) in the 
non-human primate species. Identifying the loci that 
underlie blue iris pigmentation in these species would 
help in evaluating hypotheses about the selective pres- 
sures shaping the evolution of this phenotype. 

MATERIALS AND METHODS 
Obtaining and selecting photographs 

Photographs of five black lemurs and eleven blue-eyed 
black lemurs from the Duke Lemur Center (DLC), Dur- 
ham, North Carolina, were kindly provided by David 
Haring (DLC). Photographs of 42 Japanese macaques 
housed at the Primate Research Institute (PRI) in 
Inuyama, Aichi Prefecture, Japan, were taken by PRI 
staff. We obtained photographs of 19 free-ranging Japa- 
nese macaques in the Choshikei Monkey Park on Sho- 
doshima Island, Kagawa prefecture, Japan. Photographs 
of eight brown spider monkeys were kindly provided by 
Rebecca Rimbach at the German Primate Center, 
Gottingen. Photographs of 119 humans of European 
ancestry were kindly provided by Esteban Parra and 
Melissa Edwards of the University of Toronto, Missis- 
sauga. Detailed information about photographic methods 
may be found in Supporting Information Section 1. We 
removed any photographs that had low resolution, were 
too dark to recognize iris color by eye after white bal- 
ance adjustment, or did not include the complete iris, 
leading us to exclude nine captive Japanese macaques 
and two brown spider monkeys. For individuals with 
multiple photographs of sufficient quality, we randomly 
selected one for analysis. 

Processing and summarizing iris pixels 

We applied white balance to each photograph using 
Adobe® Photoshop® automatic white balance adjustment 
and selected the iris with Adobe® Photoshop® CS4 quick 
selection tool. We excluded all pixels with RGB values 
higher than 250 (black) or lower than 10 (white), which 
represent shadows cast by the eyelids and eyelashes or 
reflections of light in the iris. We used the CIE L*a*b* 
color space, in which L* represents brightness, a* repre- 
sents relative amount of magenta to green, and b* repre- 
sents relative amount of yellow to blue, for further 
analysis (Malacara, 2002). We converted from RGB to 
CIE L*a*b* using custom code by Mark Ruzon down- 
loaded from http://www.mathworks.com/matlabcentral/ 
fileexchange/24009 on January 25, 2012, and calculated 
the median of a* and b* (excluding L* because of its de- 
pendence upon lighting conditions). Further description 
of the choice of white balance method, consistency of iris 
selection, and choice of color system and summary 
statistic may be found in Supporting Information 
Sections 2 to 4. 

Test of two clusters versus one 

To test a two cluster model of the phenotype distribu- 
tion against a one cluster model within each primate 
group, we fitted a mixture normal model with two com- 
ponents and compared the likelihood to that of a simple 
normal model using a likelihood ratio test (Supporting 



Information Table SI). Models were fitted using the mix- 
tools version 0.4.4 (Benaglia et al., 2009) and fitdistrplus 
(Delignette-Muller et al., 2010) packages in R 2.13.1, 
and an empirical null distribution of the test statistic 
was generated by calculating test statistics for 100,000 
datasets simulated from the simple normal model using 
the maximum likelihood estimate (MLE) of parameters 
(see Supporting Information Section 5). In order to mini- 
mize the effects of relatedness in the lemurs, we 
included only minimally related individuals (no relation- 
ships closer than avuncular), resulting in a sample size 
of eight. In order to reproduce the effects of this limited 
sample size in humans and Japanese macaques, we 
selected eight individuals by randomly sampling the 
minimally related individuals, and we performed tests of 
clustering on 1,000 such random samples. 

DNA samples 

Blood or DNA samples were obtained for eight of the 
photographed Japanese macaques housed at the PRI. 
DNA was extracted when necessary using the QIAamp 
DNA mini kit (Qiagen). Fecal samples were collected non- 
invasively from nine of the photographed Japanese maca- 
ques from Choshikei Monkey Park. Initial visual categori- 
zation of these macaques suggested that four (adult 
females) had blue irises, one (adult female) had irises of 
intermediate color, and four (two adult females and two 
juvenile males, one of which was the offspring of the adult 
female with intermediate iris color) had brown irises. 
Samples were stored in ethanol and transferred to silica 
gel after 5 to 12 h. DNA was extracted using the QIAamp 
DNA Stool kit (Qiagen), using the standard protocol with- 
out stool tubes. Four samples (892, 893, 902, and 903) 
were incubated in buffer ASL for 10 min at 70° C using a 
heat block in step two to increase DNA yield. 

Blood was obtained for six black lemurs and six blue- 
eyed black lemurs from the Duke Lemur Center (DLC) 
during annual check-ups or before euthanasia. Two 
black lemurs (Blanche-Niege and Louie) and two blue- 
eyed black lemurs (Bogart and Lamour) were also 
included in Bradley et al. (2009). DNA was extracted 
using the QIAamp DNA mini kit (Qiagen) and whole ge- 
nome amplified using the GenomiPhi DNA amplification 
kit (GE Healthcare). 

PCR and sequencing 

Primers were designed to amplify subsets of the con- 
served region orthologous to the region surrounding 
SNP rsl2913832 in humans (Supporting Information 
Fig. SI and Table S2), using the available sequences of 
Lemur catta (GenBank AC 126425) and Macaca mulatta 
(GenBank NC_007864). Conditions for PCR protocols are 
indicated in Supporting Information Table S3. PCR 
products were purified using Exo-SAP. Sequencing reac- 
tions were run using BigDye® on an Applied Biosystems 
3730XL or 3130 capillary sequencer (see Supporting In- 
formation Section 6). 

Alignment 

Sequences were aligned using eBioX vl.6bl (http:// 
www.ebioinformatics.org) with the default parameters of 
the Kalign algorithm. Alignments were adjusted man- 
ually when this created a consensus among all sequences 
by eliminating gaps. A M. mulatta sequence (GenBank 
NC_007864) was used to scaffold the alignments from 
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Japanese macaque fecal DNA, and H. sapiens sequences 
(GenBank NC_000015.9 and AC_000147.1) were used to 
determine position relative to rsl2913832. Each base pair 
of homozygous sequence reported for wild Japanese maca- 
ques was supported by at least seven sequencing reads to 
limit the probability of allelic dropout (i.e., the failure of 
one of two alternate alleles to amplify) to <0.01 (Navidi 
et al., 1992; Taberlet et al., 1996). Sequences were consid- 
ered heterozygous when at least one clean chromatogram 
peak supported only the alternate allele or at least two 
were heterozygous (as determined by eye). All sequences 
are available on GenBank (KC693559-KC693587). 

Association tests 

We tested for association of CIE a* and b* with age, 
sex, origin (captive, i.e., from the PRI, or wild, i.e., from 
Choshikei Monkey Park) and genotype for Japanese mac- 
aques. We first tested for association with adjusted age 
(two times age divided by maximum age), sex, and origin 
separately, and we included only significant variables 
from these analyses as covariates in the genotype associa- 
tion test. To account for relatedness among individuals 
when known, we fitted linear mixed models using 
GEMMA v0.91 (Zhou and Stephens, 2012). We used the 
median of CIE a* and CIE b*, quantile-normalized to a 
standard normal distribution using R software (R Devel- 
opment Core Team, 2011), as phenotypes. The effects of 
any covariates were modeled as fixed effects, and the 
effect of relatedness was modeled as a random effect. 

In the lemurs, we performed a test of association at all 
SNPs, combining the individuals sequenced here and the 
four additional individuals sequenced in Bradley et al. 
(2009) for the SNP identified in that study. We imple- 
mented a test for a categorical variable indicating spe- 
cies by setting the phenotype value to 2 for black lemurs 
and to 1 for blue-eyed black lemurs in a mixed model 
analysis, including a random effect of relatedness, in 
GEMMA. We used the kinship package in R (Atkinson 
and Therneau, 2012) to calculate kinship matrices for 
macaques and for lemurs. 

We did not perform tests of clustering or association 
tests in brown spider monkeys due to limited sample 
size and photograph quality and possible relatedness 
within the sample. Age, sex, and known relationship in- 
formation for all non-human primate samples are listed 
in Supporting Information Table S4. 

RESULTS 

Evaluating quantitative variation in iris color 
phenotype 

Although some subset of all four species appears to 
have blue irises, it is unclear whether this "blue" trait is 
the same in all lineages, and hence whether we expect 
the degree of similarity in its genetic basis to be inform- 
ative with regard to evolutionary constraint. Thus, we 
first assessed the similarity in the derived trait among 
species. Given the success of quantitative methods in 
describing heritable human iris color variation and in 
detecting genetic associations (Liu et al., 2010a; 
Edwards et al., 2012), we assessed quantitative variation 
across species using the CIE L*a*b* color system. 

In this color space, individuals with blue irises occupy 
a highly similar region in all species (Fig. 1). When con- 
sidering both blue and non-blue irises, however, each 
species has a slightly different overall phenotypic 



distribution. Along the magenta-green (CIE a*) axis, 
lemurs occupy the most restricted region; given that our 
sample size is limited, a* variation may be greater in 
natural populations. Despite this limitation, variation on 
the yellow-blue (CIE b*) axis appears to be greater in 
lemurs than in any of the other species. Some wild mac- 
aques have lower (i.e., more green) CIE a* values than 
any members of the other species. Considering the rela- 
tively complete representation of iris color variation in 
humans and Japanese macaques in our sample, this 
may indicate that some macaques have more green 
irises than humans. Both macaques and humans have 
continuous variation between blue and brown (Fig. 
2A,D), but the relative positions of irises characterized 
as brown and intermediate differ slightly between spe- 
cies (Fig. 1). Brown spider monkeys occupy the smallest 
range of any species, which could be due to their limited 
sample size. Additionally, this species was the least con- 
sistently categorized (Figs. 1 and 2C); we suspect that 
this is a result of the difficulty of assessing color by eye 
in photographs of animals in the dense forest canopy. 
The consistently categorized spider monkeys fall within 
a similar region of the color space to similarly catego- 
rized individuals of the other species (Fig. 1). 

In the lemurs, iris color can be roughly grouped into 
two clusters, corresponding to black lemurs and blue- 
eyed black lemurs (Fig. 2B). In order to test whether a 
two cluster model fit the phenotypic distributions within 
each group significantly better than a one cluster model, 
we performed a likelihood ratio test, comparing the like- 
lihood of a mixture of two normal distributions to that of 
a simple normal model (see Methods and Supporting In- 
formation Section 5). This test indicated a significantly 
better fit of the mixture model to the black lemur and 
blue-eyed black lemur CIE b* data in both RAW (a mini- 
mally processed image format; P = 0.0010) and JPEG (a 
more processed image format; P = 0.00049) data, sup- 
porting the presence of two distinct clusters (Supporting 
Information Table SI). For CIE a*, the mixture model 
did not fit the data significantly better than the simple 
normal model in the lemurs (P = 0.75 RAW and P = 0.11 
JPEG). In humans, the mixture model fit the data signif- 
icantly better than a simple normal model for both color 
dimensions (a* and b* P = 0.000010 JPEG and RAW); 
the two modes in humans may represent the influence of 
the main effect SNP. In Japanese macaques, the mixture 
normal model did not fit significantly better than the 
simple normal model for CIE a* or b* (P = 0.21 and 
P = 0.31, respectively), suggesting that intermediate 
phenotypes represent a larger proportion of these sam- 
ples. This may be due to non-random sampling of Japa- 
nese macaques with a preference for lighter iris color. 
When including only eight individuals, as in the lemurs, 
none of the 1,000 random samples of Japanese macaques 
and only 3 (JPEG) or 12 (RAW) of the random samples of 
humans had likelihood ratio test statistics more extreme 
than the lemurs for CIE b* (Supporting Information Ta- 
ble S4), suggesting that the extreme clustering in lemurs 
along this axis would not be observed even in a small 
sample of the other two species. 

Association of environmental variables with 
quantitative iris color phenotypes 

In Japanese macaques, CIE a* values strongly associ- 
ated with origin (wild vs. captive; P = 0.000061). This 
effect is due to lower (more green) CIE a* values among 
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Fig. 1. Quantitative phenotypic variation in brown spicier monkeys, blue-eyed black lemurs, black lemurs, Japanese macaques 
and humans; displayed using the CIE L*a*b* color system; with color representing visible color and shape representing species/ 
population. Color of points represents the results of two people (WM and SZ) independently categorizing the irises as "blue," 
"brown," or "intermediate;" "not consistent" indicates that the two people did not agree on the category. 



wild macaques than among captive macaques, a pattern 
concordant with our observation that more wild maca- 
ques possess green or blue iris color than captive maca- 
ques. This difference could be due to population 
structure. In addition, CIE b* values of Japanese maca- 
ques were significantly associated with age (P = 0.0049): 
as age increased, CIE b* values decreased (irises become 
more blue). This result is consistent with previous find- 
ings in humans (Liu et al., 2010a). Sex was not signifi- 
cantly associated with quantitative iris color in Japanese 
macaques, in contrast to the effects found in humans, 
but this could reflect lack of power. 

Identification and association testing of genetic 
variants in candidate region 

Our sequencing identified seven new variants in blue- 
eyed black lemurs and one in black lemurs within the 
1.2 kb orthologous region surrounding the site of human 
SNP rsl2913832. Neither the newly identified variants 
nor the haplotypes (in combination with the SNP identi- 
fied in Bradley et al. (2009)) segregated perfectly with 
iris color. The previously identified variant, located 1 bp 
upstream of the human transcription factor binding 
motif, is polymorphic within blue-eyed black lemurs (at 
approximately 64% frequency) but fixed for the ancestral 
allele in the sample of black lemurs sequenced to date 
(Bradley et al., 2009; Table 1). This interspecific differ- 
ence was apparent by an association test (see Methods), 
with P = 0.00084 (Table 1). Thus, it is possible that this 
variant influences iris pigmentation, but if so, another 
locus or loci must be involved as well, given that the an- 
cestral homozygous genotype at this SNP was observed 
in one of the six sequenced blue-eyed black lemurs. 
In contrast, less than 1% of humans with blue irises are 



homozygous for the ancestral allele at the main effect 
SNP (Eiberg et al., 2008; Sturm et al., 2008; Branicki 
et al., 2009). 

We identified six variants within the 700 bp region 
sequenced in captive Japanese macaques. None of these 
variants appeared to predict quantitative iris color 
(Table 1). Additionally, no variant in the 275 bp 
sequenced was present in more than one of the nine 
wild macaques. Thus, it is highly improbable that varia- 
tion in this region strongly influences iris color in Japa- 
nese macaques. 

DISCUSSION 

Our quantitative analysis of iris color variation in the 
four primate species known to have blue irises revealed 
the blue iris phenotype to be highly similar across spe- 
cies (Fig. 1). Yet our sequencing results indicated that 
no variant within the region homologous to the causal 
regulatory region in humans is solely responsible for 
blue iris pigmentation in either blue-eyed black lemurs 
or Japanese macaques. Regardless of the genetic archi- 
tecture underlying the phenotypic variation in these 
other organisms, the genotype-phenotype map clearly 
differs from that in humans. Given that the derived phe- 
notype displays marked similarity among species, this 
implies some level of flexibility in the evolution of this 
trait. 

Our comparison of the phenotypic variation across 
species suggests hypotheses about the genetic basis of 
blue irises in Japanese macaques and blue-eyed black 
lemurs. In humans, iris color varies continuously (e.g., 
Frudakis et al., 2007; Edwards et al., 2012); however, 
our perception of this spectrum of variation has led us to 
categorize iris color into discrete classes: e.g., blue, 
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Fig. 2. Independent plots for each species or group showing quantitative phenotypic variation in the CIE L*a*b* color system, 
with color representing visible color (as in Fig. 1). (A) Japanese macaques, (B) blue-eyed black lemurs and black lemurs, (C) brown 
spider monkeys, and (D) humans. 



green, and brown (Kayser et al., 2008; Liu et al., 2010a). 
The genetic basis of iris color variation in humans 
reflects both these perceived discrete categories (>80% 
of the variation in blue vs. brown can be explained by a 
single mutation) and the underlying continuous nature 
of the phenotype (several additional loci with more mod- 
est effects have been identified; see, e.g., Sulem et al., 
2007; Liu et al., 2010a). In Japanese macaques, as in 
humans, iris color variation is continuous in CIE L*a*b* 
space (Fig. 2A). This may suggest that one or a few loci 
have strong effect, producing the perceived blue/brown 
distinction, with additional genetic or environmental 
modifiers generating a continuous spectrum between 
these extremes. In contrast, phenotypic variation is dis- 
continuous between the photographed blue-eyed black 
lemurs and black lemurs; one or more fixed genetic dif- 
ferences contributing to iris color differences between 
species could explain this discontinuity. Alternatively, 
the derived trait may result from a combination of 
changes in several genes, thus allowing specific variants 
to segregate within blue-eyed black lemurs. Additional 
phenotypic data from brown spider monkey populations 



might enable predictions about the genetic architecture 
of iris color variation in this species as well. 

The observation of a variant just upstream of the 
human causal variant in the sample of blue-eyed black 
lemurs, but not the sample of black lemurs, raises the 
possibility that this variant or one linked to it influences 
iris pigmentation. Although this variant cannot produce 
blue irises on its own, it could do so in combination with 
variation at other loci. Such a complex genetic basis for 
iris pigmentation could explain the intermediate pheno- 
type observed in hybrids between blue-eyed black lemurs 
and black lemurs (Meyers et al., 1989; Rabarivola et al., 
1991). 

The similarity in derived phenotype among humans, 
Japanese macaques, blue-eyed black lemurs, and brown 
spider monkeys makes it more likely that causal sites in 
the non-human primates lie within regions containing 
variation influencing OCA2 and other known pigmenta- 
tion genes in humans. Moreover, given the deleterious 
effects of many coding mutations in OCA2 in humans 
(see http://albinismdb.med.umn.edu) and, more generally, 
the ability of regulatory variants to produce subtle and 



American Journal of Physical Anthropology 



404 



W.K. MEYEE ET AL. 



q <H <H O ^ ^ O 



00 

I 



o 

CM 

I 



o 

I 



en 

CO 

I 



I 



CO 
CM 

ICO 
I 



o 

CO 
ICO 

I 



o 

CO 
CO 

I 



E-h 00 
AoooS<<§ 
O " ^ o 
d 



O 



.. . , tO ^ 



<! CO ^ to 

Aoo^oSSc". 

j-tj o o 



H S » - 1 <t 

q o o o 



o o ^ 



O 6 



CO CO > 

AoSoo»NS 
o o 



^ H o O "* ** ^ 



Q ^H S o ^ ^ d 



A 



^ o ^ ^ o 



> CO CO to j 

aShoo«n3 
q ^ o o o 



A<oo-<< 

q <4 O ^ ^ 



AS^-|oo-*a;^ 
^ o o o ^ 



O 



ASH 0 o^»y 
" o o o ^ 



A^ooh^S 
o ^ ^ o 



^ o ^ ^ o 



^ C+H C^ C^T Ph Ph P, 



CN ^ 

00 " 

CO CD 

00 3 

^H QH 

i-H CtJ 



CD 



t» d 

B * 

„ CD 



bo a 
o * 

o 3 



8 13 



d w 

a 

0,3 
Ph Oj 

.9 .S 

CD ^ 



JO 5_ 

~3 



CO CD 

CD CD 

=■ 5. 

a" a* 

as a! 

O O 

aj as 

S B 



0) <u 



£ 33 ^ ^ 



B-* 5 * * & 

C J -d J3 

<8 3 s 1 S 

.5 o p s t* 
^ -rH +j ^ 

CD 0) CD 
mi — CD 

3 3 



T3 T3 

CD J) 



CD CD 
CD 1 

is 

cd 53 
*5 ^ ' 



CD CD 



^ >> r? 



CD 0 

> CD 

'J 3 3. 

as o- 

, — i CD 
CD 
Ph 



s § 

cr tr 

CD CD 



cs s o 

S cd as 

B .fc "S ' 

C O 

H CD CD 

■B » 
cd as 

±5 U _H 

as ^ 

-t- 3 +3 H 

m ^ s 

■o of 

> S 1 o 

C d m 

CD CD S 

H £ Oh 



rt d 
.2 o 



3 'a 



,3 3 



tissue-specific changes in phenotype (Stern, 2000; Carroll 
et al., 2001; Carroll, 2005), regulatory regions may be a 
priori more likely to contain the causal site or sites than 
coding regions. These considerations suggest that studies 
investigating genomic regions involved in human pig- 
mentation in the three non-human primates will enhance 
our knowledge of the level of similarity in the genetic ba- 
sis of iris pigmentation across species and inform our 
understanding of the evolutionary constraint on this con- 
vergent phenotype. 
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